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ABST RACT 

The work consists of tiro parts. In part ’A ! 

parametric study lias heen made on the condensation polymerization 

of ARB type monomers and on polymerization between monomers of 

type ARA and AR*B, for the case where tlie equal reactivity 

* 

hypothesis is not valid. The analysis has been done for a 
segregated CSTR* For ARB type monomer the molecular weight 
distributions and their moments have been compared for tiro 
different kinetic models, suggested before. These two models 
have been extended to the other case and a third model is 
proposed. 

In part * B * reversible hydrolytic polymerization 
of e-caprolactam to nylon 6 in a batch reactor has been 
studied. The primary purpose was to simulate the reactions 
in the polymerization process and to establish a computationally 
efficient mathematical technique for solution. The approach 
was not successful, however, for the system considered here 
due to a problem of computational unstability. Possible 
reasons for the failure and possible ways to overcome this 
problem have been discussed. 
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CHAPTER A-1 

Int roduc t ion ; 

Usually in the study of the kinetics and the molecular 
weight distribution of linear condensation polymerizations, 
the assumption is made that the reactivity of a functional 
group- is independent of the nature and length of the chains 
to which it is ■ attached ^ “ 5 \ 

However, many systems have been reported for which this 
assumption is not valid, e.g. polymerization of alkali metal 
salts of p-chloro- and p-fluorothiophenols^ \ formation of 
polyethylene terenhthalate^ ^ and polymerization of 
diisocyanates ^ and of phenol for maldehyde systems^ . 

The violation of the equal reactivity assumption for polymeri- 
zation between acid anhydrides and aromatic diamines has been 
reported , and it has been shown^ 16 ^ that the MWE in the 
condensation product of divinyl benzene with p-cresol is 
different from that predicted using this hypothesis. 

Model studies on this change of reactivity have been 
relatively fewer. Kronstadt et.al.^"*^ used two different 
rate constants in their kinetic scheme with partial success 
in matching the theoretical prediction with experimental data. 
Miller and Macosco^ have considered stepwise poly functional 
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polymerisation where functional group reactivities in the 
.same monomer change due to a substitution effect. Ozizmir and 
Odian have analysed linear step polymerization between AA and 
tC monomers and between AA,BB and CC type monomers , 
where B and C are functional groups of the same nature but 
having different reactivity. Gandhi and Babu have considered 
the change in reactivity of functional group due to asymmetry 
or substitution effect in the molecule and due to origi- 

nally present asymmetry and induced asymmetry on the same 
molecule . Kumar, Gupta and coworkers have considered 
violation of the equal reactivity hypothesis in phenol 
formaldehyde systems ^ and have proposed two kinetic 

(21 

models on general linear polycondensation of ARB type monomers 
However, all these are studies on batch reactors. 

Condensation polymerizations in CSTR's are important 

( 2i 

industrially but few researchers have worked on it. Biesentnerger 
computed the polydispersity indices in a micromixed model of 
CSTP . Based on the same hydrodynamic model, Gupta et 
obtained the MWD and its moments for the- r two proposed kinetic 
models by them earlier for batch reactors ~^ 3 ) . 

It should be noted that for the polymerization in a 
CSTR, a. microscopic segregatidn model is often closer to the 
actual case than the micromixed model particularly because of 
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the high viscosity involved. Tadmore and Biesenlr erger ^ 
found the MWD's and their moments for condensation polymeri- 
zations carried out in a segregated CSTR, for ARB type monomer 
fo 1 lowing the equal reactivity hypothesis. In this "work we take 
the two models proposed by Gupta et al*'^ - ^'^ and make a 
theoretical parametric study on the corresponding MWD's and 
moments in a segregated CSTR. We then extend these models -which 
were originally for ARB type, to ARA + AR'B type poly condensations 
and also propose a third model to make a closer approximation 
of some real cases. In both the cases the reaction is assumed 
to be irriversible which implies that the condensation product 
is withdrawn as soon as it is produced. 
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CHAPTER 1-2 
FORMULATION 

monomer : In such a reaction at any time only one 

family of molecules will be present, namely, . The general 

K ^ 

reaction is: S m + S n -- P . >irj > ^ m+n + condensation product. The 
two models proposed are the two limiting cases of real 
polymerizations : 


M odel I 


p,11 

k^/2 


p,mn 

k^: m=n, m,n=1,2, 


= 

kp/2 ; m = n, m,n 

1 


Model II 

K p,11 = k n/ 2 

%,1n = k 11 ’ 

K p,mn = k p ’ 
= k p /2 , 


n > 1 

m ^ n; m,n > 1 
m = n, m,n y 1 


( 1 ) 


( 2 ) 


The corresponding rate equations for species and moments, in 
dimensionless form, for hatch reactors, are given in Table A1 . 


In an SCSTR there is a residence time distribution and 
in an ideal CSTP, it varies from zero to infinity with a RTD 
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function E(t) dt = - exp(-t/ *t )dt (30) (3) 

The average outlet concentration of S is given by^^’ 3 ^ 
oc n 

S n * tf'V^ exp(-t/T)at (4) 

o 

To make the equations dimensionless we define x = k S. n t a,.-' 

p 1 0 


S 


Z = 
n 


n 


aiid 9 = kS 10 t 


(5) 


'10 P 

So the dimensionless outlet concentration from the segregated 


CSTR i.s miven by 

A 

£ = A_ 
n °10 


1 

6 


oc 


z n exp(-x/$)dx, n = 1 ,2, . . . (6) 


Prom the batch concentrations of one can get the outlet 
concentration by integrating equation (6). Molecular weight 
distributions may be calculated by the following equation 

A 

nS 


mass fraction of n-mer = 


n 


A 

= nZ 


n 


(7) 


InS n 

The values of the moments in the out let of the SCSTR can 
similarly obtained by using the eauation 

A 


( 8 ) 

(9) 


A 4 f 

\ = Q \ Mfc exp (-x/0 )dx; k=0,1 ,2 ... 

where ^ n k S n / S 1Q , k= 0,1,2 ... 

The molecular weight averages can then be obtained as follows 

h = V*. . f - K'K 


( 10 ) 
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ARA + BR * B type monomers ; At any time in the reaction mass, 

three types of species •will he present viz. P , the odd A-mers; 

Q n » the odd B-mers; and S , the even-mers, as identified originally 

( 1 ) 

"by Plory . Each species will have a distribution of its own. 

The moments for the distributions are defined as follows 

P = I (2n-1 )k P n ; 1 (2n-1)* Q ; M = £ <2n) k S n (1 

n ’ =l n=.i 

The reaction scheme for the three models are shown in 

Table A-2. The rate equations for various species and moments are 

shown in Tables A-3.1, A-3.2 and A-3.3 for batch reactors. The 

values at the reactor outlet will be obtained as in equations 

(6) and (8) above, the only difference is that *x' and ' ' now 

defined as x = kMt & 0=k M *£ ( 1 2 ) 

p o p o 

In Model I, we assume that the monomer-monomer reaction 
rate constant is different from those for all other reactions, 
i.e. between monomers -polymers and po lymers- polymers . In Model II 
we are assuming that the reactivity change due to the substitution 
effect is same for both the functional groups. In Model III we are 
assuming that the reactivity of one of the functional groups 
changes due to the substitution effect, but that of the other 
remains unaffected. 
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Table A-l 

Rate equations in ARB type polymerization 


Model I 


Model II 


dZ, /dx = -(R-I)zr - M Z- 
x x o 1 


dZ 1 /dX = “ RZ X 


M 


dZ 0 /dX = | Z? - M Z„ 
Z Z JL o z 


dz 2 /ax -|z* - (R-DV 2 - V 2 


, n-1 , n-1 

dZ /dX = 4 Yi 2 Z - M Z / dZ /dX = (R-l)Z Z - + j JZ Z Z 
n z m n«m on n l n-i z £ 1 m n 

m=l m=i 


n — 3^ 4/*#* 


■(R~l )2- Z ~ M Z , 
In on' 


dM 

c 

dX 


(R-l ) ~2 % 

2 *1 2 


dM 

o 

3 >T 


n = 3, 4, . . . 

,2 


R-l ~2 


M 


2 Z 1 “T" 


(R-l)M Z. 

O I 


^ _ 

dx 


(2 k “ 1 -l)Z Z 


dM^ 

dX 


(l-2 k " 1 )(R~l)Z Z + (R-1)Z 1 


+ 2 nrfery * M i M k-i 

- M, M , k = 1,2, ... 

KL O 


k-l 


^ i : ( k ~ i ): M i 

k: 


I u (kirn" M i M k-i 


I nitial conditions : at X = O: Z n =1, Z 0 = = . . . Z^ = O? 


M» = 1, k = 0, 1, 2.,.. 


Note: For definitions of the symbols see Nomenclature 
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Table A-2 


Kinetic scheme for various AHA + BR'B -polymerization 

k 

General reactions: P + Q -JES >. n >S , „ + ¥, n,m = 1,2,... 

— — n ^ n+m-1 

P + S J ? s ’ nm > P L +¥ , n,m= 1,2,... 
n m ^ n+m 9 

k 

Q + S - Q . s - nm > Q + ¥ , n,m = 1,2,... 
n m ^n+m ’ ’ 

k 

s + S S , +¥ , n,m= 1,2,... 

n m n+m 


Model I k„„ „„ = k T ; k^ n ^ = k . n = 2,3,...; k n „* = n=2,3 


"pQ , 1 1 “ I ’ pQ , In ' ~p’ 


‘'pQ ,n1~ ‘p, J 


k pQ,nm = V n ’ m= 2 ' 5 ' • 


: ps,nm = V n ’ m = 


"Qs ,nm 


k 


ss ,nm 


k , n,m = 1,2,... 

p 

k, n^m,n,m= 1,2,... 

Jr 

k /2, n =m, n,m= 1,2,... 

Jr 


Model II 


Model III 


k pQ,1m 

• «■» 

M 

M 

H 

m = 

1 , 2 , • • * 5 

k pQ,n1 “ k I» 

n= 1,2, 

k pQ,mn 

= v 

m,n 

- 2,3,.* 

• 


lr 

ps ,1m 

= k j ; 

m = 

1 ^ 2 j • » * , 

k ps ,nm = k p’ 

• 

•» 

C\J 

a 

k Qs ,1m 

= k l5 

m = 

1, 2 , • • • J 

k Qs ,nm ~ k p’ 

n= 2,3, 

k 

ss ,n m 

= V 

n * 

m, n,m = 

1 * 2 , • • * 



k /2 , n = m, n,m = 1,2,... 

IT 


k 

k 

k 

k 


pQ,1m 

- k It 

m= 1,2,...; 

k pQ,nm " k p > n - 2 > 3 > • 

ps , 1m 

= kj, 

m= 1,2,...; 

k ps,nm = V n= 2 > 3 ” 

Qs ,nm 

■ V 

n,m = 1,2,. 

• * 

'ss ,nm 

II 

>• i 

n ^ m, n,m 

= 1,2,... 


k /2 , ri = m, n,m = 1,2,... 

Jr 


V 



Table A3.1 


9 


Rate equations for Model I in ARA + BR'B polymerization 


dP 


3-,— = - k (R-DP.Q. - k P, (fd o) + M (o) ) 
dt p 11 pl Q s 

dQ 

v ^ W, 

11 pl p 


3F 1 = - k ~ (R - 1) P -. Q 1 - k - Q i < H<0) + H i 0> > 


P 


dP_ 

S. _ v 

dt ~ p 


dQ 


n-1 


T P P -kP (f^ o) + Mi o) ) n = 2,3, 
A — r n-r p n u s ' * * 


r=l 
n 


n-1 


= *5 E P r Q n-r+l + k k p H s r s n-r ~ I k p s n *4 


(o ) 


r=l r=l 

- k p S n ( ^ 0> + ^ 0>) ' 


n * — ^ ^ g 0 + * 


d f 6 o) dff< o) , x , , 

• - -W- = - k P ^ p i °i - k P K 3 q > 


dt 

dM. 


(o) 


cRT 


dH^ k) 


s(o) s(o) 1 


f(o) 


k P (R - i)p i °i + v*p ' - i k P M <r - k P 


(o) 


(H^ o) + Mq° } ) 


ar ~ = - k P (R - 1)p i % - k P ^ k) fi o 0)+ k P E xrffcni fi E <k " 1, " (1) 


M, 


dM, (k) 


1=1 


= -k p (R-i) ~k p f^ k) m^ o) + k p Tnifrrr 


dM 


(k) 


k_ 

> 

Id 


ki. __s(k-l) s (l) 


Hq j - / h; 


a #- = 2 * k P (R - 1> p i Q i + k P E x if fcnr Q 


SrCl) pj(k— 1) 


1=0 


. 1 1 , V~ ki Sfd) »(k-l) 

+ 1 k P /- TiTk-iTi" " H s M s 

1=0 

k p S s k) e s 0) - k p B i k> < S P 0) +ff Q 0) > 


Continued . . , . 



Table A3 .1 ( con tinued ) 


1 §> 


Initial condi tion s : 


at t 

b 0 : P 1 = 

P 10 > Q 1 

_ Q 10 * s x 

V Q : 

_ / S — ° 
n' n 

, n 

= 2,3, .... 

4 0) = 

f^ 1 ' = .. 

. =^ k> = 

1 p io ? 

Si! 

O 

O 

V * 

fl 

^ — 

Q 

- = 

Q 10 ; 

4 0) 

= .. 

• =^ k> 

= 0. 


o 3 


Table A3. 2 


Rate Equations for Model II 


dp 


JL 

wr = - 

R P 

dQ, 


X 

HF" = “ 

R Q 

dp 

n v* 

t 


i ^ 0) - R p i 4 o) 
- Q i 4 o) 


ff(o) 
1 H P ' 


-k P^ H^ o) 
P n s 


dQ 


n-1 

Z 1 

r=l 

n-1 


ffl o) 

s 


at 


= k (R— 1 ) Q- s + k T Q S - k (R— 1 ) Q P, - k q 
p 1 n— l p f— i r n—r p n 1 p n 


p 

(h' o) +^ o> ) 


* 


r=l 


cH' o) ‘ = , 

ar- - k P (B - 1) E i Q i - k P (R - 15 4 \ ~ % (R - 1) p i 4 


- k P 4 0) 5 


= k p (8,1) P^ - kp (R-l) 4°\ -kp (R-l) Q l4° ^ 

- Kp 4 0) 4 0) 


Continued 
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dH 


(o) 


(R-l) P x k p (R-l) P;L Q x + k p (R— 1 ) Hp° ^ Q;l 


d t p x - ' "p 1 1 ’ ~ ”p 


+ kp ^ S<°’- 1 kp S Jo) 2 _ kp (IU1) gio) Pi _ k B 
W 0> + S Q 0> J- k p (R-l) H, 


(o) 
■p‘S 


<5 f(o) 


Q 


a^ k) 


i 

k 


at kp(R-l) P lQl + k p (R-l) Pl s ■ - 

1=1 1=1 


k ; «< i) + k 


k 
£ 


TirfelTI 4 1} ^ k “ 1) - kp (R-l) Hp k) Q! -k p (R-l ) H^ o) 
(kW(o) 


" k p * T % 




(k) 


StT 




(k) 


k 

= k (R-l)P Q + k (R-l) Q, y - . M ff(l) 

p i l . p ^ ' i z__ TTTk-ITT m s 

k 1=1 

+ k P El Trfenr S s V ^ k_1, - k p <»-i) P X ^ k> -> 

(R-l) Q 1 f£ o) - k p f^ o) M^ k) 


anr 


k k 

-2 k p (R-l) P.Q, + k p (R-l) P X ■£ TT |i_ ry , ^CU+ k 

1=0 


p 

( 1 ) 


Eo F ^ 1> «Q k-1) + k p (R-i) q x ;r trrfcn i *4 

k - 

+ * k p i rrfcnr s i 1} s i k_x> - •§ k p fi s k) s i 0> - y*- 1 * 


4 k) (p i + Q i> - k P 4 k) ' s p o) + ^ o) > 

Initial conditions : Same as in Table A.3.1. 



Table 3 
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Rate Equations for Model III 


3r = - BP i^-' p i kJ 0) 

dQ, f v 

jrr = - r p i Ql - k P Q i (p 2 + p 3 + —> - k P °i *r 


n-l 


<3p 

= RP 1 S v , 1 +k Y1 P S - k P e! o) - k P„ li o) 

at 1 n-l p r n-r p n u p n s 


dQ 


n-l 


c*T - k p Z Q r s n-r “ k P ^ Q n P 1 “ k p Q n *P " k p °n t ! 

n=l , 

dfl ( k ) - - 

air- = - k P < R - 1> p i s q o) - k P ff p k) ^ o)+ k P t*- 1 * p i 


»<o) 

S 


kl 

rm^m ivi s 


Hi 1) + 


iK *( k ) 


k ]T kj 
p fzi rmcirr 




dR Q ^ «(k) ^Ck) 

kj a(D 5(k-l) 


3 TE — - - k p n Q k) - k P ff Q k) - k p S o° ^ 0) + k p ^ 


nrV-L/ fnrV J 

xnCTrr M s 


dM, 


(k) 


dt 


= k 


<R - 1>p i £ im^hr + k P £ Tvfery; 


1=0 

= (1) g(k-l) . 1 


k— 1 

Z 

1=1 


^'^' + i k p Z rnCTn- ^s 


(X) g(k-l) 


- 3c (Erl) 


P 1 5Tg k) - k p Hg k) (g£ o) +tf£ o) ) 


Initial conditions: Same as in Table A3.1 # 
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CHAPTER A-3 


MATHEMATICAL COMPUTATIONS 

In both APB type and ARA + ER’B type polycondensations 

(31 ) 

we have used the four point Funge Kutta Gill method for 

predicting the batch concentrations. Erom these batch concen- 
trations we have calculated the SCSTR outlet concentrations, 

(32) 

using equations (6) and (8), by Simpson's rule , with a 

step size of 0.05 X . For R = 1 .0 when an analytical solution 

(33) 

was possible, a check was made by using Laguerre quadrature 
Each of the step by step integration starts at X = 0 and the 
value of each of the integrals (Equations (6) & (8)) is checked 
after every 400 step increments in X. The integration is stopped 
when it is found that even after 400 such steps the maximum change 
in the fastest growing integral is less than 0.5% of ifs present 
value. Eor ARB_ type the convergence is decided on the basis the 
integral of the second moment as this is the fastest growing 
integral. For ARA + BR'B type polymerization the convergence 
criterion is based on sum of the integrals for 2nd. moments of all 
three species. A check regarding accuracy was made by decreasing 
the step size which did not produce any significant difference. 

After stopping the step by step integration, the actual 
values of the integrals are obtained by extrapolation. This is 
possible because the rate of growing of the values of each of the 
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integrals gradually decreases as '2' increases until it becomes 
almost zero. This means that the integral values increase 
asymptotically and the final value is obtainable by extrapolation. 

The integrals for monomer concentrations gives the average 
conversion in the outlet of SCSTR . Similarly average value of the 
first three moment integrals give the number average and 'weight 
average chain lengths from which the polydispersity index can 
be calculated. 
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CHAPTER A -4 

RESULTS AND DISCUSSION 

The parameter"' used in the study of ARB type polymeri- 
zation is the ratio of the rate constants , 'R* and in the study 
of AHA + AR’B , are both ’ R' and the initial molar fraction of 

the' limiting reactant P, ’fPo'. The -weight distribution of thr- 

chain lengths for ARB polymerisation at various values of 0 

and the reactivity ratio are shown for models I and II in 

figures 1 and 2. For the first kinetic model, it is observed 

that for R> 1 , the MUD curves split into two as 'n* decreases: 

one for od^ values of n and the other, having higher weight 

fractions, for even values of n. This is expected because of 

the preferential reaction of monomer with itself thus giving 

a relatively larger rate of formation of the dimer. The un- 

roacted monomer is very high in this case. In the second 

kinetic scheme, on the other hand, a splitting is obtained 

for R<1. This can possibly be explained by the consideration 

that the reaction rate constant in dimer-dimer reaction is much 

higher than, say for exam-ole, between monomer and dimer 

(by which reaction is trimer is formed). These findings are 

(21 22 ) 

same as obtained earlier for batch reactors ’ , and 

HCSTRs^ 27,28 ^ before. 
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The number average chain lengths and polydispersity indices 
for the two models of ARB polymerization in SCSTR are plotted in 
figs 3 to 6 as function of 6 . Por the sake of comparison, 
corresponding plots for batch and CSTR are included. Also, values 
obtained with fPO = 0.5 for A HA + BR'B polymerization in SCSTR 
are included for comparison. The values for the latter ones are 
found to be always less than the corresponding ones in ARB 
polymerization. Usually the curves for SCSTR falls, as expected 
between the curves for batch and HCSTR . 

In figures 7 and 8 are compared the values of polydispersity 

indices for the three different models of ARA + BR’B polymerization 

when the initial reactant ratio is other than 1:1.’fP0' is the 

original fraction of the limiting reactant, from these plots it is 

seen that (i) the model II is most sensitive to the change in R 

values while model III is the least sensitive; this is shown by 

the fact that the change of 1 P ' is maximum for a change in R values 
to 10.0, 

from R = 0.1 : in case of model I. (ii) that In general for such 
a polymerization, the PD I value obtained for R«0 a* 1 © much higher 
compared to those obtained for R^>1 irrespective of the model and 
the fPO value. This possibly indicates that according to these 
models, for unequal initial concentrations of monomers, the 
dispersity of the polymer formed will be higher if there is a 
preferential polymer-polymer reaction as compared to the monomer- 
polymer reaction. 
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In figures 9 and 10 cross plots of J* vs monomer conversion 
for ARB type monomers, -with constant chain length curves, are 
shown. These curves are of significant interest since they show 
the relationship between the three important parameters, conversion, 
chain length and polydispersity index. Here also it is seen that 
for same conversion f falls as 'R' decreases. 

In conclusion it may be observed that though the kinetic 
models that are proposed here are idealized they may be expected 
to simulate the polymerization reaction reasonably well under 
some conditions and for many classes of polymers. The results presen- 
ted here are the outcome of a theoretical investigation and are 
given in the form of dimensionalized plots which can be used for 
study of particular classes of polymerization reactions. At least 
they can be used as a starting step while analyzing for more 
complex aspects of nonequal reactivity, provided more data are 
available . 





d> 

II 


-Weight distribution of chain lengths in ARB type polycondensation. 



Model I 



Weight distribution of chain lengths in ARB type polycondensation. 
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Fig, 8 - Polydispersity index vs. reactor holding time in 
ARA ♦ BRB polymerization. 
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Monomer conversion (7e) 

Rg. 10 -Cross plots of polydispersity index vs. monomer conv, including constant chain 
length curves. 
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CHAPTER B-l 


Introduction to Nylon-6 


Nylon— 6, which is the generic term for a synthetic linear 
polyamide in which the repeating unit of the chemical structure 
contains six carbon atoms, — CDNH— ( CH2 ) 5 — / is produced by 
the ring-scission polymerization of the cyclic monomer 
€ -caprolactam. The chemistry of Nylon— 6 polymerization has been 
comprehensively reviewed by Reimschueseel^^ , Polycaproamide 
can be produced from caprolactam by one of three entirely 
different routes viz,, hydrolytic polymerization, catinoic poly- 
merization and anionic polymerization. The cationic and 

/ Q /■ V 

anionic cases are only of theoretical interest' * , but the 

hydrolytic polymerization is by far the most important commercial 

process for Nylon-6 production. The actual plant practice may 

vary depending on local conditions and demands, but the basic 

( 1 7 *- 8 ) 

process remains the same * , Most of the commercial processes 

consist of a feed section where carefully controlled compositions 
are obtained containing molten caprolactam, water, acids or bases 
and other necessary additives. The bulk of the reaction takes 
place in the polymerization reactor, the product from which is 
treated for removal of unreacted monomer and low molecular weight 
cyclic oligomers either by hot water extraction or by vacuum 
evaporation. The monomer recovered is recycled and the polymer 
is sent to the processing section. 


Many variations in the process conditions and the reactor 
designs have been investigated^ 7 9 * 1 °)^ Analysis and simulation 


■(. K A NPUf 


of different polymerization reactor systems h^e^l^en, csr$; 




668J.3 * 
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out by different workers^* 9-12) based on these information 
some optimization studies have also been made^ 0 *'*’ 3 '^^. 

Besides the existing processes reported in the 

literature^'^ 0 ^ some novel operations have been patented 

( 24 -* 29 } 

recently , Many process modifications have been suggested 

also to increase the process eff iciency ^ \ to improve 
the product properties' ~ and for better removal of monomer 
and low cyclics^ 2Ci " 23 ^. 

In our present study, all the major established reactions 
involved with the hydrolytic polymerization of caprolactam 
have been taken into consideration. Based on the available 
information on reaction kinetics ard polymerization mechanism 
and partly with a semi— empirical approach/ equations for three 
sets of molecular weight moments have been developed, correspon- 
ding to the three different species present viz,, the bifunctional 
chains, the monofunctional chains and the closed macro rings. 

From these three sets of moments it should be possible to get 
all the necessary information about the monomer conversion, 
molecular weight distribution and their averages for all kinds 
of species at every stage of polymerization. 



CHAPTER B— 2 


Hydrolytic Polymerization of Nylon-6 

Kinetics and Equilibrium Studies : The course of reaction 

during polymerization of caprolactam can be followed by 

the end group determination as suggested by WiHoth^ 30 ^* 

Studies have been made on the polymerization kinetics by 

Hermans et.al.^ 3 ^, Kruissink et,al.^ 32 ^, Heikens et.al,^ 33 ^, 

Willoth (34 ' 35) , Dorr and Willoth (36) ‘, Tai et.al. i37 ' 53 > for 

hydrolytic polymerization with or without the presence of 

acids or amines. Studies on polymerization equilibrium 

and effect of process conditions on the MWD of equilibrium 

( 38 ) ( 39 ) 

polymerization have been reported by Fukumoto , Majury , 

40 ( 41 ) 

Tbbolsky and Eisenberg , Reimschuessel , Reimschussel and 
Dege (42) and Smith (43) . 

Mechanism of Polymerization : On the basis of the kinetic and 

equilibrium investigations made so far it has been established 
that the hydrolytic polymerization involves three major 
equilibrium reactions, namely, 

i) ring opening reaction, where the monomer, the caprolactam 
ring, is opened up by water molecules under proper reaction 
conditions, to form amino caproic acid, the shortest 
bifunctional chain, 

ii) poly condensation reaction, where the amino and carboxylic 
acid nd groups in linear bifunctional chains react with 

one another to form larger chains via amide linkages, 
with water as the condensation product; and 


iii) poly addition reaction, where the monomer caprolactam 

adds directly to any of the growing linear chains as in 

44 

a typical chain growth polymerization . Poly addition 
is the primary reaction during the polymerization process 
and is mainly responsible for the conversion of 
caprolactam. But in the later stages of polymerization 
linear poly condensation becomes important since this 
reaction determines the final molecular weight distribution 
in the product, as equilibrium is approached. 

Besides these three, two more reversible reactions 

have been identified. One, the so called termination, is 

the reaction of bifunctional linear chains with monofun ctional 

acids or amines, which are most commonly used as chain 
(45 V 

stabilizers , to form monofunctional long chains. The other 
reaction, known as cyclization, is a transamidation by which 
the end group of long poly caproamide chains attacks its own 
amide linkages, by what is termed as 'back biting', to form 
shorter linear chains and macro rings of moderate molecular 
weights. This reaction will again be discussed later; 

Cyclic Oligomers : In the equilibrium polymer of Nylon-6, some 

low molecular weight macrocyclic compounds are always found. 
Though their concentrations are rather low, the removal of these 
compounds from the polymer is essential prior to processing, 
since their presence adversely affects the processing of the 
polymer and its physical properties. For this reason these 
compounds have been subject to considerable investigations. 
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detection 

Their / and quantitative estimation has been done by 
many workers^ 47 ' 53 ^ . Mochizuki and Ito^ 11 ^ have reported 
cyclic oligomer concentration as a function of reaction time 
during polymerization. Based on the ring- chain equilibrium 
concept of Jacobson and Stockmayer^ 4 \ Semiyen and Walker^^ 
had attempted, using the Flory-Williams rotational isomeric 
state model, to predict theoretically the percentage of low 
cyclic oligomers in Nylon-6 equilibrium melt, though the 

i 

comparison with experimental data is inconclusive, Andrews 
et al^ 4 ®^ have reported experimentally determined equilibrium 
concentrations of cyclic oligomers up to the hexamer and the 
corresponding equilibrium constants for cyclization reaction. 

They also described the effect of post polymerization treatment 
on their concentration, Reims chuess el ^ has attempted to 
curve fit these equilibrium constants as a function of the 
ring lengths and temperature. 

Rate and Equilibrium constants ; Reports on the investigation 
of reaction kinetics of caprolactam polymerization have shown 
that each of the ring opening, poly condensation and poly addition 
reactions are catalysed by -CDOH end groups ^ 31 ' 33 \ though 
the actual dependence of rates on catalyst concentration is not 
established rationally, in the usual empirical approach the 
forward reaction rate constants of these reactions are thought, 
as a first approximation, to consist of two parts, one of Which 
is linearly dependent on -ODOH group concentration, and the other 


independent as follows* 

k i * + c k* 


( 31 ) 


( 1 ) 



where 1 i* designates either of the above three reactions, 

*c* is the acid end group concentration; k? and k^ are similar 
to Arrhenious reaction rate constants, each being determined 
by a frequency factor and an activation energy e.g. •' 

k ± = A^ exp (- e| / RT) (2) 

where superscript * j * is either * o* or 1 c* . ihe equilibrium 


constants K*s for these three reactions have been determined 


from thermodynamic quantities 'Ah' and 1 4S' and many reports 


agree quite well. The relationship is 

As. AH. 

K i * exp ( — (3) 


lb completely specify the kinetics, one thus needs to know 
A i* E i # ^ S i 311(3 aH i as 9 iven in Reimschuessel^ 1 ^. Recently, 
however, Tai et al ' have determined these values independently^ 
using least squares curve-fitting and though their reported 
reported numerical values sometimes differ substantially from 
Reimschuessel 1 s values, they are claimed to produce less 

i 

deviation from their extensive experimental data. 

Mochizuki and in their model, have tried to 

show that the forward reaction rate constant in the reaction with 
monofunctional acid depends on -CDOH concentration as well as 
the melt visoosity, though their data show that the viscosity 
dependence is quite small. In our case, we will consider that 
the reactivity of an end group in a monofunctional chain, for 
example, is the same as that in a bifunctional chain. 

In the cyclization reaction also, Mochizuki and I to ^ 3-3 ^ 
have taken an entirely different approach which is not quite 
unquestionable. They have used, for a cyclization reaction. 


o / 


the reaction between the -NH 2 and group and -ODOH end group 
of the same chain, with production of one water molecule and 
a macroring of the same lengths 

e.a- S C_ + n = 2 , 3 ,,,. (4) 

^ n-*= — n z 

They did not consider the attack on mid chain amide linkages 

by end groups, which is the general cyclization reaction* 

Moreover, their data for the rate constants are not published 

due to proprietory reasons* In the present case we have used 

a different cyclization reaction scheme and as a matter of 

fact, we have not considered the above reaction of Mochizuki and 

I to at all* The reason is if we consider the above reaction 

with the same theoretical approach as in the other reactions, 

it is imperative that we consider the ring opening of 

macrorings also, as it is a reversible reaction. But, except 

Tai et«al„^ 3 ^ who considered ring opening of the dimer, no 

data are available in the literature for such a reaction.' 

Hence, we have kept it out of consideration. However, we can 

partially justify it by considering the fact that the 

probability for this reaction would always be the least among all 

other cyclization reactions for the same chain length because it 

involves complete end to end bending of the polymer chain. In 

our model we have been forced to make certain simplifications on 

(54) 

the original model of Jacobson and s tockmayer in order to 

make it amenable to mathematical treatment. In a semi— empirical 
approach we have taken the cyclization equilibrium constant to 
vary inversely with the length of the macroring formed. We 
will consider it in some detail in tfie next chapter* 


MWD : 


Theoretical studies on the molecular weight distribution 

(57 ) 

in poly caprolactam has been relatively fewer. Willoth and 
( 58 ) 

Hermans et.al. have tried to match the MWD in poly caprolactam 

to a Flory-S chul2 distribution with relatively little success, 

this is more so in the later stages of polymerization. Tirrel 
(6 ) 

et.al. simulated the polymerization of caprolactam, in 
different PFR-CSTR combinations, under some simplifying assum- 
ptions (and with some errors) and reported the theoretically 
predicted MWD averages. Tai et.al.^ 59 ^ recently tried to 
predict MWD averages theoretically with an assumption of 
S chulz— Zimm distribution, and found polydispersity indices s 

from GPC about 9 % higher than computed valves. 

In our work we will be making use of an approach 
suggested by Krishnan^ 12) and Gupta and Gandhi (61) which is based 
on the method of Min^ 60 ^. Krishnan has simulated some 
CSTR-PFR combinations for Nylon-6 polymerization and ha s reported 
the MWD averages, while Gupta and Gandhi have suggested a general 
method for getting the actual MWD as well. The mathematical method 
used by them is computationally efficient, but in none of the 
above works the cyclization reaction has been included, which we 
will be doing here. 

In this context, references to two other recent works 
may be made, though they are not particularly related to Nylon-6. 

In the first, Howell et. al.^ 62 ^ attempts to get a semi empirical 
distribution function for linear and circular Gaussian chains; 
and in the second, Siootte^ 63 ^ uses a Gram-Charlier series 
expansion of molecular weight distribution. This latter approach 
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is in many ways similar to that used by us, but comparatively 
more difficult computationally and there is no evidence that 
it is guaranteed to work for Nylon-6 system. 
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CHAPTER B— 3 
Reaction Simulation 


Reaction Scheme : The present reaction scheme which includes 

all five reactions identified so far in Nylon-6 polymerization 
was first proposed by Gupta et.al.^ 6 ^. The reaction scheme is 
presented in Table B-l. At any time in the reaction mixture, 
there would be three kinds of species; e.g . bifunctional chains, 
denoted by S , n=l,2,...? monofunctional chains, denoted by 
A^, n = 1 , 2 ,... , and rings, denoted by C , n = 1 , 2 ,.... 

Initial feed is aomposed of * C^' (i.e. the monomer)? water, 
designated by *W', and 'A^ 1 , the stabilizer (usually monofunctional 
acid, e.g. ROD OH, where 1 R* is a nonreactive group). The rate 
constants k 1 , k 2 , k 3 , k£, k£ and k^ are based on the concen- 
tration of the functional groups like -000H, -nh 2 and — CDNH— 
etc. and these are the ones reported in the literature^ ^ . 

In writing rate equations for molecular species, we have to 
account forall the reactions oc curing. 


i) 

ii) 


for ring opening reaction, rate equations are directly 
obtained for molecular species, e.g. c^. 

constant 

for condensation reaction, in the forward reaction/would 

be 2)^, as there are two possibilities, the -CDOH of 

may react with -NH 2 of S m and vice versa, unless n and m 

are equal, when it would be k 2 ^ 64#65 \ in the reverse 

reaction, the rate constant will usually be 2k^, because 

a water molecule can attack a linear chain S , at two 

n-fm 
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places to give exactly similar products? the rate 
constant will however be k£ if (n+m) is even and the 
central amide linkage is attacked, 

iii) for reaction with monofunctional acids, the forward 

reaction rate constant is always k 2 , as only one 

-NH 0 on S may react with one -COOH of A . The reverse 

z n m 

reaction rate constant will always be kj,, because the 

reaction between W and A gives different products 

n+m 

when water attacks the i-th and the (n+m-i)th amide lin- 
kages, 

iv) in the general polyaddition step, since the - NH^ end 
groups of S n can react with any of the ' m* amide 
linkages in a cyclic m-mer, the rate constant is mk^. 

In the cyclization step, the — NH_ end group of S ■ 

* " u i> 

can attack any of the (n+m-1) -OONH- groups on itself 
to give several products? the rate constant for this 
reaction depending on the size of the macro-cyclic 
formed (C^) and denoted by k^ (i). 

The cyclization reaction rate constants require 
special consideration here. First of an let us establish 


the relation between this cyclization reaction and that 


originally dealt with by Jacobson and Stockmayer 


( 54 ) 


i.e . 


n+m 


C + S 
— n m 


( 1 ) 


It diould be noted that these two reactions are of opposite 
nature? i,e. the way the reactions have been defined, the 
forward reaction of the first is equivalent to the reverse 
reaction in the second, the difference being that the former 
reaction considers equilibrium between all the chain sand rings at 



a time, while the latter concentrates on a single chain in 


equilibrium with a ring of a particular size* It has been 
shown'- that the equilibrium constant for reaction (1), 
K a , defined as 

K = ( s *c\/s ( 

n V m n / / m+n 


can be expressed as 


K n = V*" 


(3) 


where *x* is the fraction of reacted end groups in the chain 
fraction. Now if we consider our present cyclization equilibrium 
it is easy to establish that 


oo 


n 


nk. 


m— 


n+m 


jfm if — »■»*■■■ 


m 


k 3 ( n) ^ 


m 


m=l 


1 /, Xitel 

= CL * V “ CD 

n e 

m=l 


-> 


m 


(4) 


Defining an equilibrium constant Kgfn) for the present case 
as Kg(n) = k 3 /k^(n), we get 


K 3 (n > - 


( 1 — 


r s ra 

m=l 


n 


oo 


n 


oo 

E 

m=l 


.) = 


m 


nK^ 


(1 


E sVZ 

m=l rf _ , 


m 


m=l 
(5) 


It should be noted that the concentrations and 

*x' are the equilibrium values. Hence x is very nearly unity. 
So from eqn. (5) it may be concluded that at least for small 
values of ‘n*, it may be written, wihout much error, as 

K 3< n > - ttV* <6) 

Now for n s 1, is the equilibrium constant for the reaction 
+ S n s n+ i* ihis value has been investigated by many 
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workers and agree reasonably well ^ ^ . 

The values of K for n = 2 to 6 are reported in 
n 

the literature ^ . values for n y 6 are not available because 
it‘s not possible to detect and quantitatively measure the 
concentration of in equilibrium polymer beyond n y 6. 

From the data of K n vs. n , it is shown that fortunately the 
functionality can often been expressed as 

/ n * 2,3,.... (7) 

where B isa constant. As an example *B' is approximately 
0.1692 at 235°C (64) . Thus from (6) and (7), one gets 

k * B 

(n ) = * n s= 2, 3, . . . (8) 

We have assumed that this functionality between k^(n) and n 
is retained for values of n greater than 6. As will be 
seen later, this functionality is a very convenient form which 
enables us to adopt the moment approach. For any functionality 
relationship other than that in (8), a moment analysis becomes 
extremely difficult. 

Rate Equations ; The reaction mechanism having been established, 
we may now proceed to write the rate equations for various 
species present in the reaction mass. The rate equations for 
a batch reactor or PFR are shown in Table B-2. For a CS TR, 
the derivatives on the LHS of these equations will be replaced 
by the difference in input and output concentrations divided by 
the reactor holding time. Because of the reversibility of 
the reactions* it may be seen that the solution for any species* 



say s - will require knowledge of the concentrations of- 

n 

all three species A,c and S, longer than 'n'. For this reason, 

the analysis becomes, at least theoretically, the solution of 

an infinite set of equations. Were it possible to solve all 

of these coupled differential equations (or algebraic equations, 

as the case may be), all the necessary information about the 

/ \ 

polymerization could have been obtained. Gupta et.al. 
have integrated this set of equations for batch reactors and 
have reported the results in the form of molecular weight 
distributions and their averages, for both chains and rings, 
as function of reaction time. However, as is expected, this 
approach takes disconcertingly a large computation time, which 
is quite natural because of the necessity of handling of literally 
hundreds of coupled differential equations. This has motivated 
search for more computationally efficient mathematical approaches 
even at the expense of some accuracy, if necessaryl The moment 
approach which will now be described is undertaken with this 
view in mind. This approach will presumably be more useful, 
in case of CSTR's, since it will be extremely difficult, to 
solve so many algebraic equations efficiently. 

Moment equation ; -As has just been stated, to make a rigorous 
and accurate analysis of the polymerization reaction, one needs 
to handle, at least theoretically, an infinite number of 
equations. The idea in the moment approach is to convert this 
infinite set of equations to a set of reasonably small number 
of moment equations so that the analysis becomes easier. 
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As has been previously described there are three kinds 
of species in the polymerizing mixture* bifunctional chains (S ' 
mono functional chains (A^ ) and rings (C n ); and each of these 
species will have a weight distribution of its own. 

Let us define three sets of moments as follows : 

- E nk v\ = J= 

n=l n “ 1 

CO 

- E (9) 

n=l 

If we multiply the rate equation for any species, say 
S , by n , and sum up the resulting equations for n=l to oo , 
we should get the rate equations for d ^^/dt. The same 
method applies for the ^‘s and C n * s* After considerable 
mathematical manipulations one arrives at the set of equ ations 
given in Table B-3. Note that no expression is available 
for d _/dt, because that involves summation of a harmonic 

° a*, 

series. Also note that — - gg - + — ^ ^ = 0 , which is 

ejected since + #1 + ^ ]_ = constant = + a Q * 

•To integrate the rate equations for the moments it is 
necessary to get explicit rate equations instead of those 
involving series summations as shown in the above equation. 

In order to do it, recurrence formula for these series summation:- 
are necessary. Appendix B-l lists some of these series sums. 

Samples of individual moment rate equations are shown in 
Appendix B~2. Using the series summation formula listed in 
Appendix B-l, any other moment rate equation can be similarly 
obtained. 


A, 


and 

l 1 
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It should be noted that the rate equation for the k-th 
moment requires the knowledge of the (k+l)th moment. Our 
subsequent approach is to search for a mathematical tool using 
Which one can express the (k+l)th moment in terms of the k-th 
moment, and get the actual molecular weight distribution from 
theoe moments. This method will be described in the next 
chapter. 

Once we are in a position to express the (k+l)th 
moment in terms of the k-th moment, we can integrate the above 
set of differential equations until equilibrium is reached. The 
results thus generated will give all the necessary information; 
and after a sound mathematical method is established, one may 
possibly use this technique for an efficient reactor simulation 
purpose and for subsequent optimization studies. However, 
that is beyond our scope. 



47 


Table B-l 

Kinetic scheme for nylon 6 polymerisation 


1 . Ring opening : 



s 


1 


2 . Polycondensation : 

2k 2 <n7$n) or k 2 (n=tfn) 

forward reaction : S n + S s n+m + 

2k 2 

reverse reaction : S n+m + W > s r + s n+rri ~r' 

r — 1/2,... 

Note : The above rate constant will 

be ki if (n+m) is even and 
2 2. 

for r = ^-( n+m ) 
ink ^ 

3. General polyaddi tion : + C m 9* s n+m • n ' m * 

kl (r) 

Cyclization : S n+m > S n+m«r + C r' 

r — 1 , 2 , . . « 

Note: The above two reactions are 
also valid between A-species 
and C -species 


4. Reaction with mono- 
functional acid : 


S + A 
n m 


-s*" A + W , 
n-H-ft 


n,m = 1 / 2 /... 


Reverse 


A + W 
n+m 


k 2 


"9* S Hh A ✓ 

r n+rn—r 

2C 22 1 2 / # * 


Contd # * * 



Table B-l (continued) 


Structures : 


n 


H- 


H 


-N - 


(ch 2 ) 5 


°7 

«7 


C 


► OH 


n 


A n : 


R- 


T~C 

it 

C 


H 

I 

N 




(ch 2 ) 5 - 


0 

II 

c 


- OH 


n-1 


C : 
n 



C 1 %s the 


monomer. . 
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Table B— 2 

Rate equations In isothermal nylon 6 reactor 

(for brevity, concentration terms are not enclosed by square 
brackets ) 


ds^ 

dt* 


k-jC^W - k|S 1 - 2k 2 S 

oc 


«« OC 

1 , s n + 2 *2 W 2- s n - k^ Z nc n 

n=l n=2 n=l 




Z. k«(n) s_, - k 0 s n Z A n + k£W 


n=l 


n+1 2 1 


n=l 


n=2 


n 


ds 




n 


n-1 


dt 


- 2k^S X- S + k„ 21 S S ■ - k'W(n-l)s 
2 n - m 2 -*~r m n-m 2 n 

m=x m=i 

^ n-1 

+ 2k;,W 2. S - k S 2. mc m + k 3 21 <n-”> s m C n-m 

m=n+l m=l m=i 


n-1 




s_ ^ k ^ m) + X k*(m) s_ m - k 0 S 
m=l 6 m=l 




n -5—? -3 N “ W T r~; ''3' 11W ~n+m ~2 n m 


dC-j 

d t 


+ k 2 W ^ . *ta 
m=n+l 


- kjC w + kp i 

+ k’(l) 5" A 
J n=2 n 


n = 2, 3, 4 


ex: ^ 

CL s + ki(l) 2- S - k„c, 

i 1 n n 3 „ « n 3 1 , n 

n=l n =2 n=i 


dC 

n 

St, 


k~nC 2 4 - k* (n) S 


— 

3 n -1 ffi 
m=i 

+ k^(n) Z- A m 
m=n+l 


m=n+l m 


k 3 nC n *Ti ** 
m=i 


Contd 
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Table B 


dt~ = 


d \i 

anr = 


■2 ( continued ) 





<C*C 






I 

n=l 

S 

n 

+ kiW 

z n=2 

A - k 
n 

A 2 
m=l 

IllC + 4^. 

10 rnl 

m=l 

k‘ (m )A 

3 n+rr 


€> <Z 


1 






- k 2*h 

z 

m=l 

s 

m 

11 

+ k 2 S \i S n-m 
m=l 

- k^W(n- 

- 1 ,A n + 

k 2 

W A 

— . ir. 
m=n+l 


<=*c 


«"C 






- k 3^ 

E 

m=l 

mC + k_ <2l_ 
m 3 ~r 

m=i 

(n-m)S, 

C 

m n-m 

n m=l 

k 3 

(m) 


oc 

X 

m=l 


k' (m)A 
3 rn-m 




dt 


V" - 44 + k 3 V^ "Vlb ' b) Qo - ,s i > 


aA k 

■at" 


-\V-* 2 \p o + T ^2 ^ 

K 

k MW - k'S - 2k % % + k 2 n , (k k ' i n ' t T 

m=o 

— k^W( + + 2 + »•♦.+ (m—1 ) ^ 


+ (m-1 ) v S 


“ k 3 ^k'^1 + k 3 V^ "ITW “ li t ^ k~l^ : 


k.(i- - B> 2. W* 1 - k(k “ 1} m k “ 2 + ... + i-U M )S r 


*X C“C 


. _ V '!> k-l . -k(k-l) „k-2 _ . . , ^k~l^ 

-•k^B / - n- - -\kfti- " — 2 — m n 4- . . .+ t-n J > s m 

n=l m=n+l 


k 2?>k(V +k 2 W 


2 *fl k % 2 k + ... + (m~l) 


4 MM + 4*1 " k 3/\3k + l + k 3 < £T " B)< ^0 “ s i 


Oc o^c 


+ k 0 B -. 


n=l m=n+l 


2. 2, V k4 + 2 k “ 1 + ... + <m>-l ) k->1 / S 


Cont<3 



trz 



+ k 3 ( K. 


oc 

B)((? - 2L) +k,B 2 

^ n=l m=n+l 


’l + 2 krl + 


-f (m-1) 


A 


m 


<3/3 


VO 

"eft" 


' k 2 + k 2 fo^o " k 2 W( fl ~ (V + k 2 W( (*1 ~(V 



k 3 Po ’ 1 + k 3 P o 'l - ° 


“at" 


■ k 2 + k 2 +E SirfesTT^m 'V-m " k 2 w( pk+l ”(V 


+ k£W 


<5 ( 

i { 

m-9 


m=2 
k 


l k + 2 k + ... + (m-l) k £ ~ k,/^^ 


k 3 ll(k-i)l fk-l~^ 


m 3 
zxr 


i+i - k 3<lr - B) ^ (kmk_1 


k 3 m=2 

o > c 


k(k-l ) m k™2 , . / i .. t- -n ^ 5 (km k ”’’*' 


2 m + ...+ (*-!) )A m “ k^B 


n=l m=n+l 


k(k-l) k-2 
_i_ m n + . 


. + (» n ) k *" 1 )A 


m 



CHAPTER B-4 


Computational Procedure 

We have seen in the last chapter that the solution for 
the Nth moment requires the knowledge of the (N+l)th moment. 

Our next task is to (i) find the (N+l)th moment from the lower 
moments and (ii) to interpret the N moments in order to get 
the actual mo lecular weight distribution. We will aJopt Min's 
approach^ 60 ^ as has been done before by Krishnan^ 12 ^ and Gupta 
and Gaidhi (61) . 

The starting point in this method is the consideration 
that any continuous normalized distribution G(x ) can be expre- 
ssed by the aim of an infinite number of Laguerre polynomials 
as follows; 

co 

.jf e -x v a. « 

G(x) “ p W+T) ~ 0 C n L n (1) 

where L n (x) is the generalised Laguerre polynomial expressed 
as 

hi (x) “ H ( - 1)j “31 (2) 

j=0 

and ^ is a constant parameter. 

The reason for using Laguerre polynomials is its more 
general orthogonal properties and the capacity of taking care 
of fractional moments as suggested by Min^ 60 ^. Though theoreti- 
cally an infinite number of polynomials are necessary for all 
practical purposes; it is possible to express G(x) by a reason- 
ably small number of polynomials. Here we have taken the upper 

limit of summation to be N w = 9, The actual choice of N _ 

max max 

depends on the storage capacity of the computer, the accuracy 
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desired and the computation speed wanted* Once a finite number 


is used, because of the orthogonal properties, the unknown 
«r 

coefficients C n , n = O to N max' 0311 ^ computed if the first 
N max moments of the normalized distribution G(x) are available 


a 

"n 


n 


L&+11 nl y / 1} J 
r(iHoi+i) -C- 1 


j=o 


r (n -Wf 1) G j 

fTj +*+l)(n-J)i Ji 


(3) 


Thus if the first moments of G(x) are known* 

<K 

C n * s can be calculated therefrom and the actual 
distribution, G(x), can be expressed with reasonable accuracy. 
Once this is done, the (N max + l)th and highermoments of the 
normalized distribution can be obtained using the general 


expression 


N 


na+k+i) 

r<4(+ i) 


r(ot+ i) 

P( n -K(4- 1) P(k+j +6(+ 1) 
(n - j)i ji fTj+^+ 1) 


max 

X 

n=l 


o( 

r* 

^n 


i <-i)i 


j=0 


(4) 


for derivation of (3) and (4) see Appendix (B.3) 

The basic approach is thus to convert a given set of 
say *k' original, unnormalized moments to a set of ' k* 
normalized moments. Next by the method described above, the 
higher normalized moments are obtained, from Which is 
regenerated the (k+l)th moment of the original distribution. 

lb convert the original moments to a set of normalized 
moments, the first task is to simulate the actual discrete 
distribution by a continuous distribution. Let us designate such 
a distribution by M(v)dv, which has to be converted to normalized 
distribution G(x) dx« Obviously it follows that 
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M(v) dv = M^ o) G(x) dx 
oo 


(5) 


since J[ G(x)dx is unity, 
o 

Min has suggested that there be a linear relation 
between v and x e.g, 

v = a.x. (6) 

the constant * a* has to be determined from other conditions. 
Since *&* is a free parameter in eq.(3), it is usually customary 

Q< 

to put Cj as zero* In such a case it follows from eq.(3) 


that 


G 1 = 1 +<X 


Now from eg* (5) we can write 


OD 

S 


v* 


M (v ) dv = M. 


(o) 


GO 

I 


or 


u v 


M^ o) a k §< k) 

V X 


(7) 


v G (x ) dx for any value of 


( 8 ) 


Putting k = 1 in eq. (8), we get 

* 

v 

a = — 7 y V = - 

h(o) 


M. 


Cl) 


V 


and 


(l+«) M, v 


M (k) = 

* V y 


V 


(1+C0 E^ o) 


(l+oO M^ o) 


pr(k) 


(9) 


(10) 


Eqn.(10)relates the original moments to the normalized 
moments. It should be noted that this eqn. works only for k=2 
and above. Por kpl eqn. (10) simply becomes eq. (7). It only 
shows that *(X* and G^ are in ter- related but neither has a 


unique value. 

u 
C. 


To fix a unique value of 'OC'/ Win has suggested tha 
2 be also put to zero. For such a case, eqn. (3) yields 

j 

G 2 b (1+(X) (2 4-00 B G x (1 + G x ) (11) 



56 


the other hand putting k = 2 in eq. (10) yields 




( 12 ) 


Combining eqs. (11) and (12) one gets 


G 1 - 


h (o 

u y V 


(13) 


Thus G^ and hence get uniquely '.determined. 

Once the G(x) distribution is obtained, M(v) distribution can 
also be obtained by making use of equations (5) and (9)j yielding 


M(v) ss 


H^ o) (l+CX. ) 


°< x > fi< o) (H0C) 

X = 5<1) " V 

% 


(14) 


lhus for integral values of v, m(v) can be calculated? 

and though, strictly speaking, M(v) is a density distribution 

function and not the actual distrifc* tion, i.e. actual discrete 

b* bar 

distribution of say S n is related as s n = J* M(v) dv* 

n^a. a 

to calculate individual concentration S n , the evaluation of 
M(v) at that value of v is sufficient. Because Av = 1 in the 


integral, and (M(v) can be takei to be equal to S n# 
remarks apply to all other distributions* 


The same 



CHAPTER B— 5 


Computational Problems and Discussions 

There are two main computational problems present 
in this case. Firstly, it was seen that though the value 
of a suggested by Min gives faster computation, it does 
not necessarily give good convergence. The convergence can 
be best tested in the following way. While using Min*s 
procedure we start initially from a set of N number of 

ITtctAi 

unnormalized moments from which G^'s are calculated. Now 

we may recalculate the normalized moments by putting k = 2,3, #tt 

♦ „„ in Eq.(4) in Chapter B-4, which, when unnormalized, 

max 

can be compared to the original set. The sum of the deviations 

from the original moments constitutes the best convergence 

t 63} 

criterion. The same has been suggested by Sicotte'' also. 

For a particular set of original moments, this sum of 
deviation is a function of cc only, and it was seen that in 
most of the cases, the a suggested by Min gives a poor 
convergence. The value of a that gives the best convergence 
may be obtained by an optimization process, where the sum of 
deviation is to foe minimized. This a is subsequently used to 
compute the (N maX +l )th moment. 

The second problem comes while applying Min 1 s method to 
the ring length distribution. It should be noted that there 
is no way to directly find the zeroth ring moment, since 
it necessitates the summing of a harmonic progression, for 



which there is no convenient formula. But the value of 
zeroth moment is essential in order to apply Min*s method. 

Hence one has to think of some way to circumvent these 
problems. 

In the present case we have used a as a truely free 
parameter for the chain distribution and obtained a by an 
optimization as justified above, one surprising finding is that 
integer values of a always give better convergence. This 
practically makes it an integer optimization process. However/ 
it is felt that more studies are warranted in this direction. 

In case of rings/ a has again been used as afree 
parameter. is obtained from Eq.(7) and the zeroth moment 

calculated from G ± by Eq.(13)/ followed by the optimization 
process some as in case of chains. The difference in the 
procedures for chains and rings is that in the former G^ 
is calculated from the zeroth moment and in the latter. It 
is just the opposite way. This is possible presumably 
because there are two inherent degrees of freedom in this 
computational scheme. 

Using the approach outlined above we have tried to 
solve the thirty coupled differential equations which arise 
in the reaction simulation viz. the rate equations for the 
zeroth and first nine moments of both the chains species 
and rate equation for monomer and the rate equations for the 
first to ninth ring moments using a four-point Runge-Kutta-GIll 
method ^ 66) with an initial time increment of 0.1 hour. The 
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rate and equilibrium constants of Reinschuessel^^ are used 
for an initial water concentration of 0.44 mole/kg and 
temperature 235°C* 

The molecular weight distribution for the chains can 
be found by the method outlined in chapter four and they 
match within some errors with the data given by Gupta et al.^^. 
The same method to find the ring length distribution for the 
macrocyclics was not very successful, however, particularly 
at the low values of the polymerisation time. This is 
presumably because the monomer concentration is so much 
more than the total concentration of macrocyclics that the 
distribution is not very well defined (it should be noted that 
in the calculation of ring moments the monomer concentration 
is also included). This problem may be circumvented by a 
proper scaling procedure in which the monomer concentration 
is deducted from each of the ring moments. This, understandably ^ 
won't instill any significant error because the contribution 
by the monomer concentration is the same to all of the ring 
moments. By this scaling process a reasonable distribution 
may be obtained for which the concentration for n—1 should be 
ideally zero. In fact, minimization of the concentration of 
n=l after scaling in itself may be selecting criteria for a , 
but in this case the least squared deviation may be poor. It 
should be admitted that we are yet to establish a really 
sound selection criterion. 



But the most difficult problem which could not bo 
overcome in the present case and which was the cause of the 
failure . was the stability problem. That is to say. the values 
of dYj^/dt for larger values of k become so highly negative 
after about two hours of polymerization, that Y ^ ' s themselves 
gradually became negative and the computation becomes 
unstable. This un stability is typical of nylon-6 polymeri- 

/ z: >7 \ 

zation as has been reported earlier by Gupta ot al* » 

In their case they used individual species equations and 
commented that the problem may be overcome by sufficiently 
decreasing ^t . But in our case we could not get over the 
problem even with a value of = 0*001 hr* This also put 

into doubt the acceptability of the method in which we are 
computing the zeroth ring moment. This may well be one of 
the reasons for this instability. The zeroth moment can be 
evaluated, for example by curve fitting the harmonic 
progression by some suitable polynomial. But whatever may be 
the reason, unless we are able to solve this problem on a 
sound theoretical basis, it is not thought worthwhile to go 
on decreasing to solve the unstability problem. 

In conclusion it will be said that here our aim was 
mainly to establish a generalised technique so that the simu- 
lation con be done for any set of operating and system variable 
For this particular case, we have not been successful due to 
problems outlined above, but it is hoped that, with more in- 
depth study, .these problems can be solved. This method can be 1 
subsequently used for better reactor designs and optimization. 
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Appendix A-1 


Defini tion 


S n-1 t k = 

S n-1 ,1 = 

S n-1,2 = 

S „ , = 
n-1 ,3 

S n-1 , 4 

S n-1 , 5 = 

S , c = 
n-1 , 6 

S n-1 , 7 = 

S n-1 ,8 = 

S - n = 
n-1 ,9 


1+2 k + 3 k + .... + (tt-l) k 


n(n-1 ) 

2 

2 

2n -3n +n 

6 

„4 0 3,2 

n -2n +,n 

4 

6n 5 - 15n 4 + 10n 2 -n 
30 

2n 8 -6n 9 +5n 4 -n 2 
12 

12n^-42n 8 +42n 9 -1 4n 9 +2n 
84 

3n 8 -1 2n 7 +1 4n 6 -7n 4 +2n 2 
24 

1 0n 9 - 45n 8 + 60n 7 -42n 3 +20n 3 -3n 
90 

2n 10 -10n 9 ± 15n 8 -14n 6 + 10n 4 -3n 2 


20 


Sample rate equations for the moments 


Bifunctional 

moments 


Monofunction 
Momen ts 


chain 


d A 8 

“ 3 RT 


cs k n MW 


k 

11 


k 2 ( 16^^ 7 + 




+112X^5 + 70 X\) - 4 W -A 8 ) 

10^+45^ + 6oA ? - 42 + 20^ - 2% x 
+ 2 k 2 W ( r 1 """ 1 " 1 1,11 ■ r " r 1 ^q' " ^ 

+ h 2 ie^j' i 2 + 2Q\% + 56^5% + 70\% + S6^ 3 9 6 


Chain 


+ 28 A 2^7 + 8 A 1 ’ 9 8 + ^o ' : ?9 5 

- h 3 (— - B)( 8^ 7 ^ 28^ 6 + 56^ 5 - 7 °A 4 + 56 A. 
3 

-28> 2 + a\ + 0 - S x ) 

-_k 3 B -6(^g -? 7 ) + 56(— ^ - - - ■ ) - 70 

A" 2 VV , H _ 15^ + 1°^ 

* 5 > + 56 < 30 > 


2 q( 2 Ab - 6 ^ 7 + 5 \“\ 


12 


•)+ 8 ( 


12 ^- 42 ^ 4 - 42 ^- 14 ^+ 2^2 

, 


3> 8 - 12 A ? + 40> 6 - 7 A 4 + 2^ 2 
.( ) 


24 


1 0 ^ 9 - 45 | 3 - q + 60 ^- 42 /? 5 + 20 /? 3 - 3 ^_ 


- k 2 ^(?o + k i W( 
df B 

" ^^oVflV 28 (VV^/VV 
+ S6 h\ + 28 (W 8 p7\>- k 2 W % - P 8 > 
+k 2 W( 


l 0 p 9 - 45 ^ 8 + 60 ^ 7 - 42 ^ 5 + 20 | 3 3 - 3^ 1 

- 


-)+k, 


(8(3 7 i 2+ 28^ 3+ 56(J 5 4 4+ 70 ( 3 4 - ) ? 5+ S6,«39 6+ 28 ^ 7 
+ 8 fl ? 8 + fo 1 (’ 9 , - k 3 (_ i^ -B ) ( 8 ^ 7 - 28 ^ 6 + 56^ 5 - 70 ^ 
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a^~ - - k i MW +k i s i- k 3 ; \>V k 3 ( -^ - :: i 

+K 3 B( 24 > 


x , 3 Pq*" 12 /V +1 * ; f fe~ 7 *W2 

M “ 53 a 


- k 3| 3 oV k 3 t ’K; - B)( , 5 o - V +k 3® 
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APPENDIX B-3 


MIN * S METHOD 


G(X) 


a -x 

x o. 


N 


a 


max a 

21 V*> C „" 


n a+1> n=0* 

“ , , £ (_1)J rtn+a+l) _4. 

where L n (x) = ' ' (n-j)! f"(j+a +l) j! 


N 


How, 

or 


a G “ x oc 

, G(x) (x) = n a + 1) L n (x) ^ C m L m (x) 


niox (x a 
m 


m-0 


N 


n * /y tnt /y rnscjc 

X ft °<*> ■ V*> ^ <£ 


i=0 


m=0 


Integrate 

co 


n •? 

( ST ( i%J ,n+a, . _ F( n+« + 1) J* 

j ^ (- 1) <n-j } J r m ^ r( a +1) ni C n 

do i=o 


X=0 






OO 


NOW 


5f) 



G(x)dx = 


max X M e -x 
n p(a +lT 


n 

j=o 


,n+a, xi 

ji 


)dx 


or rearranging/ 
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